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The calculation of scalar gravitational and matter perturbations during 
\ multiple-field inflation valid to first order in slow roll is discussed. These fields 

I \ may be the coordinates of a non-trivial field manifold and hence have non- 

minimal kinetic terms. A basis for these perturbations determined by the 
^ ■ background dynamics is introduced, and the slow-roll functions are general- 

ized to the multiple-field case. Solutions for a perturbation mode in its three 
^ \ different behavioural regimes are combined, leading to an analytic expression 

5j \ for the correlator of the gravitational potential. Multiple-field effects caused 

by the coupling to the field perturbation perpendicular to the field velocity 
can even contribute at leading order. This is illustrated numerically with an 
example of a quadratic potential. (The material here is based on previous 
work by the authors presented in |hep-ph / 01 07272| . ) 
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1 Introduction 



As has been known for a long time, inflation offers a mechanism for the production of 
density perturbations, which are supposed to be the seeds for the formation of large scale 
structures in the universe. This mechanism is the magnification of microscopic quantum 
fluctuations in the scalar fields present during the inflationary epoch into macroscopic 
matter and metric perturbations. Also, since a part of the primordial spectrum of density 
perturbations is observed in the cosmic microwave background radiation (CMBR), this 
mechanism offers one of the most important ways of checking and constraining possible 
models of inflation, see e.g. [Q, especially when combined with large scale structure data 
i- 

There are two important reasons for considering inflation with multiple scalar fields. 
The problem of realizing sufficient inflation before a graceful exit from the inflationary 
era and producing the observed density perturbation spectrum in a model without very 
unnatural values of the parameters and initial conditions can be solved by the introduction 
of additional fields. This is the motivation for hybrid inflation models ^. The other 
reason is that many theories beyond the standard model of particle physics, like grand 
unification, supersymmetry or effective supergravity from string theory, contain a lot of 
scalar fields. Ultimately one would hope to be able to identify those fields that can 
act as inflatons. In addition such string-inspired supersymmetric models naturally have 
non-minimal kinetic terms. 

The previous two paragraphs outline the motivation for looking at perturbations in 
multiple-field inflation. A lot of work in this direction has already been done, for example 
in D, I, |, |ll|, |12|, ^ Q (more references can be found in [l^). However, most of the 



previous literature does not consider the most general case, but is usually limited to only 



two fields and/or minimal kinetic terms. The papers [|T^, |TT[ do treat the general case, 
but the authors do not consider the effect of rotations of the basis, nor do they work out 
explicitly the particular contribution to the gravitational potential. In |T5| we provided 
a general treatment by computing the scalar gravitational and matter perturbations to 
first order in slow roll during inflation with multiple real scalar fields that may have non- 
minimal kinetic terms. Which of these fields acts as inflaton during which part of the 
inflationary period is determined automatically in our formalism. 

This paper basically summarizes part of our previous work ||15| , concentrating on the 
calculation of the gravitational potential during inflation, in particular on multiple-field 
effects like the influence of entropy perturbations during inflation. Necessary background 
concepts, like the induced field basis and the generalized slow-roll functions, are also 
discussed. The result is presented in the form of the gravitational correlator at the time 
of recombination, but for the evolution after inflation only adiabatic perturbations are 
considered here. This paper also contains a new result for the spectral index n — 1 of the 
perturbation spectrum. 
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2 Background 



The gravitational background of the universe is described by the flat Robertson- Walker 
metric with scale factor a: 

ds^ = -dt^ + a(t)2dx2 ^ ds^ = a{T]f {-drf + d^^) . (1) 

The comoving time t and conformal time i] are related by dt = adrj. Derivatives with 
respect to t and t] are denoted by a dot and a prime, respectively; the associated Hubble 
parameters are H = a/a and Ti. = a' /a = aH. Another useful variable is the number of 
e-folds N, defined by a{N) = aoexp(A^), which can be considered as a time variable as 
well: dA^ = Hdt = Hdr]. 

For the matter content of the universe we assume an arbitrary number of real scalar 
fields, which are the components of a vector (p. Apart from a generic potential V(</)) 
we allow for the possibility of non-minimal kinetic terms, encoded by a field metric G. 
In other words, the scalar fields may be the coordinates of a real field manifold with a 
non-trivial metric G. This is a common situation in for example supergravity models, 
where the (complex) scalar fields parameterize a so-called Kahler manifold with a metric 
that is the second mixed derivative of the Kahler potential. Since a complex field can 
always be written in terms of two real fields, our treatment with real scalars is sufficiently 
general to be easily applicable to these special manifolds. 

The Einstein and field equations lead to the following equations of motion for the 
homogeneous background: 

i^' = y + H = -\k'\<I>\', A0 + 3if0 + G-iV^y = O. (2) 

Here we have defined the covariant time derivative on a vector in field space as VtA"- = 
A" + T'^^^'^A^. (Indices a,b, . . . are used for components in field space and F^^ is the 
connection associated with the field metric G.) The V is used for covariant derivatives 
with respect to the fields: {'W)a = dV/d(j)°'. The length of a vector is given by |A| = 
VA ■ A, with the inner product defined by A • B = A"''B = A^GB. The quantity n is 
the inverse reduced Planck mass: = SvrG = Sn/Mp. 

The background field dynamics induce a prefered basis on the field manifold. The first 
unit vector ei is given by the direction of the field velocity 0. The second unit vector 62 
points in the direction of that part of the field acceleration Vt4> that is perpendicular to 
the first unit vector ei. Hence: 

ei = ^, e2 = ^^^, P^ = l-pH = l-e,el, (3) 

with pll the projection on the direction (f). This Gram-Schmidt orthogonalization proce- 
dure can be continued to construct the remaining basis vectors [^, but the first two are 
the only ones we need here. 
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We can now define the following multiple-field slow-roll functions: 



,^.JL f,^3± i^m. 

of which we can take components with respect to the basis defined above, for example 
r]" = ei ■ f), f/-*" = e2 ■ r) and ^" = ei • j^. Again, one can easily define slow-roll functions 



of higher order [15|. With these definitions as short-hand notation we can rewrite the 



background equations of motion in the following form, which is still exact: 



K ^ / 1 



The assumption that the slow-roll functions e and ff are (much) smaller than unity is called 
the slow-roll approximation. (The second order slow-roll function ^ is assumed to be of 
an order comparable to e^, er^H, etc.) If this assumption is valid, we can use expansions in 
powers of these slow-roll functions to estimate the relevance of various terms in a given 
expression. For example, to first order the Friedmann equation (^ is approximated by 
replacing (1 — e/3)^^/^ by (1 + e/6). The background field equation up to and including 
first order is given by (^) with the right-hand side set to zero, as all those terms are order 
3/2 or higher. 



3 Perturbations 

On top of the homogeneous background treated in the previous section there are small 
quantum fluctuations. We consider only scalar perturbations and write the matter and 
metric perturbations as follows [0 

0f"^r/,x) = (/>(r^) + (50(r/,x), 

C\V, x) = 9,M + S9,AV, x) = a'ir^) (^"J - 2a\vMv, x) . (7) 

All equations are linearized with respect to the perturbation quantities. The gravita- 
tional potential $(17, x) describes the scalar metric perturbations and is a quantity we 
are interested in, as it is related to the temperature fluctuations in the cosmic microwave 
background radiation by means of the Sachs- Wolfe effect. 

Instead of using 6(f) and $ to describe the perturbations, it turns out to be more 
convenient to use the so-called generalized Mukhanov- Sasaki variables q and u, as well as 
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a short-hand notation 6, defined by 



u 



These redefinitions remove first order time derivatives, so that it is easier to understand 
the behaviour of the solutions, and are necessary for quantization. The perturbation 
equations in terms of spatial Fourier modes k now read ||T5| 



(6"\ 
k^-jUk = nfj^e^ ■ qk, 2^?qk + {k^ + ^'f^)qk = 0, 

f2 = -^M^ - (2 - e) U - 2e (^(3 + e)P'l + e^f,^ + f/ej ) . (9) 

Here the effective mass matrix is defined by = G'^V^VV - R{(j), 0), with R the 
curvature tensor on the field manifold, [R(B, C)D]'^ = R^^^B^C^V^. From its definition 
we can derive that 9"/9 is exactly given by e"/e = n'^{2e + t)" + 2{fi'^\f - {fj^f - i}). 
We see that in the multiple-field case the redefined gravitational potential u is coupled 
to the field perturbation in the e2 direction. In the hterature (see e.g. [|]) perturbations 
in the ei direction are called adiabatic perturbations, while perturbations in the other 
directions are called entropy perturbations. This coupling of the gravitational potential 
to the entropy perturbations is suppressed by the slow-roll function r^-*-, so one could a 
priori expect the contribution of this inhomogeneous term to be only important at first 
order. However, because of integration interval effects it turns out that it can contribute 
even at leading order, as shown in the example in section ^. 

To quantize the perturbations, we explicitly choose to work in the basis {ei, e2, . . . } 
defined in the previous section, denoting vectors in that basis with non-bold symbols: 
1^ = (in) = (©n ■ q)- Although this has several advantages, most importantly it results 
in a standard canonical normalization of \{q')'^q' in the Lagrangean, independent of the 
field metric G, making quantization straightforward. Of course there is the price that 
terms with Z (see below) appear, but these can be dealt with. We can write qirf) = 
Q{ri)a} + Q*{7])a, with constant creation and annihilation operator vectors and a and 
a matrix function Q{ri) that satisfies the classical equation of motion. Finally we perform 
a rotation Q = RQ to simplify this equation of motion and find: 

Q" + ik^ + n^n)Q = 0, R' + nZR = 0, = 1 ■ I),e„, (10) 

ri 

with Cl = R~^QR. The matrix Z is anti-symmetric, first order in slow roll and only non- 
zero for m = n±l. The presence of this matrix is caused by the fact that the basis vectors 
are in general not static in field space because of the background dynamics. Working out 
the definitions, one can easily show that Z21 = fj-^; general expressions are given in |p!5[ . 

Considering the equations and (11^) for u and Q and reahzing that TC grows rapidly 
during infiation, while k is a constant per mode, we see that their solutions change dra- 
matically around the time rj-^ when a scale crosses the Hubble scale ('passes through the 



4 



horizon'), defined by T-C{rin) = k. Hence there are three regions of interest: sub-horizon 
(7Y <^ A;), transition {Ti ~ k) and super-horizon {Ti. ^ k). (Notice that we are considering 
a single, though arbitrary, mode k here, on which the resuhing expressions will depend.) 
The equations in the sub-horizon region are easily solved, but the oscillatory solutions 
there are irrelevant for the correlators we are interested in. 

In the transition region we can only solve the equations analytically if we make the 
additional assumption that the slow-roll functions are constant. Since the derivatives of 
the slow-roll functions are one order higher in slow roll, this is a consistent approximation 
to first order, provided that the transition region is small enough. The solution for Q 
valid to first order in slow roll in the neighbourhood of tj-h, taking into account the correct 
initial conditions, can be written in terms of a Hankel function: 

Q{v) = ^Hjil{kv), un = lt + 5n, Sn = int - ^ + 2~en e,el (11) 

Here the additional assumption has been made that also those components of the matrix 
M^/H^ that cannot be expressed in terms of the slow-roll functions defined in (^, are 
nonetheless of first order in slow roll. In particular this puts constraints on the curvature 
tensor of the field manifold. The matrix Z appearing in equation (10) for R is a possible 
source of multiple-field effects. However, it turns out that in generic situations the effects 
of the rotation of the basis during the transition region encoded by Z are beyond first order 
so that it does not appear in the first order result (pJ^, although it might be important if 
fj^ peaks around r^T^. 

Solutions for u and Q can also be determined analytically to first order in the super- 
horizon region. However, matching the solutions in the super-horizon region to those in 
the transition region to determine the constants of integration is not trivial, as there is 
no region of overlap where both solutions are valid to first order. Hence the standard 
method of continuously differentiable matching at a certain time scale cannot be applied. 
In particular one should not simply match at rj-j-c, as the assumption of neglecting k, that 
is necessary to obtain the analytic super-horizon solution, is not valid there. 

We solved this problem in by setting up series expansions in kr] for both the 
super-horizon and transition solutions. We found that the leading order powers in these 
asymptotic expansions can exactly be identified to first order, for both the decaying and 
non-decaying independent solutions, which can already at zeroth order be distinguished 
from each other. Since the transition and super-horizon solutions are approximations of 
the complete solution of the same equation of motion, the coefficients in front of the series 
have to be equal as well. In this way the solution in the super-horizon region is determined 
completely to first order. 

The final result for u at later times in the super-horizon region (i.e. neglecting terms 
that are suppressed by the scale factor) to first order in slow roll is ^ 



n dn' n dn' H' 

uM=Up^ + D^e -L, up^ = 9 -L dr]"n9fj^q,^, (12) 
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with Dk = — I ^.J^ dn^i ((1 — ^h)^ + B6n) o) + c.c. and up the particular solution 
caused by the inhomogeneous source term in the equation of motion for u and 
5 = 2 — 7^; — In 2^ 0.7297. Switching back to the real gravitational potential $ we 
have the complete solution at the end of inflation, so including the effect of the perpen- 
dicular field or entropy perturbations during inflation. Considering here for simplicity 
only adiabatic perturbations after inflation we can compute the vacuum correlator of the 
gravitational potential at the time of recombination when the CMBR was formed ||15||: 



|2 _ //is2\ _ 9 



+ IB {{2ln + 4) + '^nn^l^Pe + UIM^p^ ■ (13) 
Here Up^ = Up{rie), with r/g the end time of inflation, is defined by 

[/J = 2v^ r dv'n^^ ^ e^gg^i (14) 

and we have used the fact that Up has no component in the Ci direction. More information 
on Up, in particular on how to rewrite it in terms of background quantities only, can be 
found in . Apart from this amplitude we can also compute the slope of the spectrum. 



This spectral index can be determined one order further, but most important is its leading 
(first) order part: 

mn\6£_ . 2Ul^{2~en + fli-5n)Up.-Af,^elUp^ 

^-^hT^- ^''^^ l + UlJJp, • ^^^^ 

The explicit multiple field terms in the results for |(5kp and n — 1 are the contributions 
of the terms with Up^ and r/^, which are absent in the single-field case. Since they are 
both to a large extent determined by r^-*-, we see that the behaviour of f/-*- during the last 
60 e- folds of inflation is crucial to determine whether multiple- field effects are important. 



4 Example 

In this section the example of a quadratic potential on a flat field manifold is briefly 
discussed to illustrate the theory. In this case it is possible to compute Upe explicitly 
and we show that it can indeed contribute already at zeroth order to the gravitational 
correlator, despite the slow-roll factor fj-^. 

The potential is \^ = cji^m^cf) with a general symmetric mass matrix given in 
units of the Planck mass . With a trivial field metric and initial condition </>(0) = 0q 
we can write the general first order slow-roll solution for cf) in terms of a single scalar 
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Figure 1: a) The slow-roll functions and b) the particular contribution Up to the gravitational 
correlator as a function of the number of e-folds in the model with two fields on a flat 



manifold with a quadratic potential with masses rrii 
conditions (pi = 20, 02 = 25. 



1 ■ 10- 



m2 



2.5 ■ 10"^ and initial 







Contribution 


Error 


n — 1 


Contribution 


Homogeneous 


1.55 • 10"'^ 


0.505 


0.0001 


-0.0377 


0.584 


Particular 


1.52 • IQ-^ 


0.495 


0.0006 


-0.0269 


0.416 


Total 


3.08 • 10-y 


1 


0.0003 


-0.0646 


1 



Table 1: The amplitude |5kp and slope n — 1 of the vacuum correlator of the gravitational 
potential are separated into a purely homogeneous and a (mixed) particular part. Given are 
their values and their relative contributions to the total according to our analytical slow-roll 
results (p^ ) and (|T5|) combined with (0). For the amplitude we also calculated exact numerical 
results and give here the relative error between those and our analyical approximations. 



function V'(^)- Even without knowing this function exphcitly we can determine Upe HTSj : 



m = e-i-^^W0o, t/Je = (16) 

This expression for Upe assumes slow roll, but is valid at the end of inflation provided 
that fj^ goes to zero. 

For the case of two fields with masses mi = 1 ■ 10~^ and m2 = 2.5 ■ 10~^ and initial 
conditions (pi = 20 and (p2 = 25 in Planck units, the results for the slow-roll functions 
and Up as a. function of the number of e-folds are plotted in figure |I|. We see that in 
this case fj-^ indeed goes to zero at the end of inflation after reaching a maximum during 
the last 60 e-folds, and Up reaches a constant value long before a possible break-down of 
slow roll. Values for the gravitational correlator amplitude |(5kP and spectral index n — 1 
are given in table |I|. 

The results for the amplitude and the slope are split into a homogeneous part (all 
terms without Upe) and a particular part (the rest, including mixing terms). Everything 
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is evaluated for the mode k that crosses the horizon 60 e-folds before the end of inflation. 
From the column giving the relative errors between our first order analytical results ( ]T3| ) 
and (|16D on the one hand, and the exact numerical result on the other, we see that 
these results agree with our claim that we computed the correlator to first order in slow 
roll: the relative errors are (much) smaller than (9(e-^). We also see that our slow-roll 
approximation for Up is indeed still very good at the end of inflation, as indicated by 
the small error in the particular part. From the other columns we see that the particular 
solution terms are responsible for almost half the total result in this model, both for the 
amplitude and the spectral index. Hence neglecting these terms to leading order, which 
might naively be done because they couple with a fj-^ in (P), can be dangerous. 

5 Conclusions 

In this paper we have given a summary of our general treatment for scalar pertur- 
bations on a flat Robertson- Walker spacetime in the presence of an arbitrary number of 
scalar fields that take values on a curved field manifold during slow-roll inflation. These 
are the kind of systems that one typically obtains from (string-inspired) high-energy mod- 
els. Here we concentrated on the calculation of the vacuum correlator of the gravitational 
potential to first order in slow roll, which is related to the temperature fluctuations that 
are observed in the CMBR. 

A discussion of the background served as the foundation for this analysis. The back- 
ground field dynamics naturally induce an orthonormal basis {ei, e2, . . . } on the tangent 
bundle of the field manifold. This makes a separation between effectively single-field 
and truly multiple-field contributions possible and is a necessary ingredient for correct 
quantization of the perturbations. We also modified the definitions of the well-known 
slow-roll parameters to define slow-roll functions in terms of derivatives of the Hubble 
parameter and the background field velocity for the case of multiple scalar field inflation. 
These slow-roll functions are vectors, which can be decomposed in the basis induced by 
the field dynamics. For example, the slow-roll function fj-^ measures the size of the field 
acceleration perpendicular to the field velocity. Because we did not make the assumption 
that slow roll is valid in the definition of the slow-roll functions, it is often possible to 
identify these slow-roll functions in exact equations of motion and make decisions about 
neglecting some of the terms. 

We generalized the combined system of gravitational and matter perturbations of 
Mukhanov et al. |jl6| by defining the Mukhanov-Sasaki variables as a vector on the scalar 
field manifold. The gravitational potential only couples to the scalar field perturbation 
in the direction e2 with a slow-roll factor fj-^. First order solutions were found in the 
three different regimes that reflect the change of behaviour for a given mode when it 
crosses the Hubble scale. Using a procedure of analytically identifying leading order 
terms in asymptotic expansions in ki] we were able to relate the solutions in the different 
regions and find the complete first order result for the gravitational potential at the end 
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of inflation. Considering only adiabatic perturbations after inflation we could give this 
result in terms of the vacuum correlator of the gravitational potential at the time of 
recombination. We also determined the spectral index n — 1. 

Multiple-field effects are important in the adiabatic perturbations at the time of recom- 
bination if fj-^ is sizable during the last 60 e-folds of inflation. The most important source 
of multiple-field effects is the particular solution of the gravitational potential caused by 
the perpendicular field perturbations during inflation. We found in our numerical example 
of multiple scalar fields on a flat manifold with a quadratic potential that this term can 
contribute already at leading order, even though it enters with a slow-roll factor in the 
equation of motion. This is true for both the amplitude and the slope of the spectrum. 
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